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Abstract 

We introduce a dynamical matrix model where the matrix is interpreted as a Hamiltonian rep- 
resenting interaction of a bosonic system with a single fermion. We show how a system of second- 
quantized fermions influences the ground state of the whole system by producing a gap between the 
highest eigenvalue of the occupied single-fermion states and the lowest eigenvalue of the unoccupied 
single-fermion states. We describe the development of the gap in both, strong and weak coupling 
regime, while for the intermediate coupling strength we expect formation of homolumo "kinks". 
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In the complex systems composed of fermions interacting with bosons, such as a molecules, 
nuclei etc., appearance of the energy gap between highest energy level occupied by fermion 
and lowest unoccupied level, so called homolumo gap [if, is a well known fact confirmed 
by experiments as well as by exact calculations from the first principles in various specific 
examples. Some of the most important properties of the system, stability, interaction with 
another system, size etc., are determined by the physics in the neighbourhood of the gap. 
In a more general setting, we are interested in the application of the homolumo effect in the 
project of Random Dynamics 0]. There, one starts from the observation that the energies at 
our disposal are extremely low, compared to the fundamental energy scale, presumably to be 
identified with the Planck scale. Consequently, from the fundamental scale point of view, the 
usual high-energy physics can be described as low-energy excitations in the neighbourhood 
of the gap. 

The purpose of this letter is to provide a model which possesses the main features of 
the interactions involved in the production of the homolumo gap and allows generalisations 
which incorporate other properties such as the appearance of the single level in the gap or 
mixing of the densities of the occupied and unoccupied levels. This letter can be viewed as 
continuation of our previous work |3| which, using different approach, gives a better insight 
into mechanism of production of the gap. In this new approach we would like to confirm that 
the homolumo gap arises whenever we have a system of fermions and bosons in interaction, 
provided that bosons are sufficiently soft to yield to the pressure from fermions, while the 
details of the model itself seems to be unimportant. 

The main assumption of the model is dynamics of single fermion energy levels. This 
assumption can arise, for example, from the observation that the energies of the single 
electron levels change as the nucleus of a molecule vibrates. We show that the Hamiltonian 
describing the system of interacting bosons and fermions can be constructed from very 
general assumptions and can be written^ as 

H = Hb + 9Hfb, (1) 

with Hb and HpB given by 

Hb = uj{B^B + 55t), Hfb = pBf + /t^t/. (2) 

This form of the Hamiltonian appeared in description of black holes j3] , mesons and hadrons 
in QCD [sj] and Jahn- Teller effect [l|. Within our approach, the part of the Hamiltonian 
which describes fermion-boson interaction is completely determined, while the Hamiltonian 
describing bosonic degrees of freedom is a matter of choice. Our choice is motivated by 
possible applications in different branches of physics. As written, the Hamiltonian Hb has 
been used in description of (a sector within) supersymmetric Yang-Mills theory in four 
dimensions, two-dimensional quantum gravity and string theory \d\, and matrix cosmology 
0]. It is also related to integrable models such as the Calogero model Q . Furthermore, the 
ground state wave function of Hb corresponds to the Gaussian ensemble from the Random 
Matrix Theory (RMT) Q which has been successfully applied in analyses of the spectra of 



complex molecules and nuclei, transport properties of disordered mesoscopic systems |10 



^ The equations should be understood as matrix equations; B's are represented by matrices and /'s as rows 
and cohimns. 
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We start modelling the Hamiltonian by considering system of the interacting fermions de- 
scribed by the Hamiltonian in which we include appropriate antisymmetrization of fermionic 
modes in analogy with symmetrization of bosonic modes: 

Hfb = \Y. {fif - ffl) Xl (3) 

l,k 

where matrix X is, for the moment, fixed NxN hermitean matrix and //, z = 1, . . . , are 
fermionic operators satisfying the usual anticommutation relations. Using the appropriate 
unitary matrix U, the matrix X and the fermionic operators can be rewritten as 

xi = Yl ul'^^ui, fi = j2 dlut. f = Yl ul'd\ (4) 

i k I 

where Xj's are eigenvalues of the matrix X. Expressed in the new coordinates the Hamilto- 
nian HpB becomes 

HFB=lj2^^{did'-d'diy (5) 

i 

As the Hamiltonian Hps does not change the number of fermions, the eigenstates can be 
chosen with definite number of fermions 

|^^n,....J=<---4|0)/, (6) 

where |0)/ represents fermionic vacuum with property rf^lO)/ = 0. As usual, we interpret 
operator dl as creation operator for a fermion at level i with associated energy Xj and 
|^ji,...,i„) describes state in which levels ii, . . . ,in are occupied while the rest of the levels 
are unoccupied. Therefore, the form of Hps is completely determined by the requirement 
that the eigenvalues Xi are interpreted as the energy levels of single fermion. 

Next, we introduce a model where the matrix elements Xj are not fixed, but are random 
variables, distribution of which is determined by a probability law P(Xj). Then, we can 
naturally ask the question about densities of occupied and unoccupied levels in a certain 
state of the system. These are given by the expectation values of the operators 

N N 

p°-(x) =J2S{x-Xi)dld\ p"(x) =J2S{x-Xi)d'dl. (7) 

i=l i=l 

For the fixed matrix X, the densities ([7]) in a generic eigenstate are given by the appropriate 
sums of the delta functions. Particularly, in the ground state of the system, the interval of 
non-vanishing density of occupied levels lies below the interval of non-vanishing density of 
unoccupied levels. The main consequence of introducing the distribution law P{Xj) is that 
the sharp delta-functions profiles are smeared, thus allowing the penetration of the density 
of the unoccupied levels into the interval in which density of occupied levels is nonzero and 
vice versa. In such situation one expects that the gap between highest occupied level and 
lowest unoccupied level disappears. We explore this question in the setting where matrix 
X is dynamical matrix, so that probability law P{Xj) is a consequence of the quantisation 
of the dynamical degrees of freedom^ of the matrix X. In that case the matrix degrees of 



^ One might say that we approximate fundamentaUy random physics by more tractable Hamiltonian dy- 
namics. The other point of view would be to attribute the success of RMT to the underlying Hamiltonian 
dynamics. 
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freedom contribute to the total energy of the system and to be specific, we assume that 
self-energy of the matrix X is determined by the following Hamiltonian: 



HB=ij:pF^+Y^^Fi^ (8) 

where Pj = —id/dXj is momentum conjugate to Xj. After quantisation, Hb can be ex- 
pressed in terms of the appropriate bosonic operators and together with HpB, which we 
now interpret as Hamiltonian of fermion-boson interaction, constitute the system governed 
by the Hamiltonian ([1]). The parameter g defines fermion-boson interaction strength and 
can be taken as positive number without loose of generality, since for negative g we can 
transform X — >■ —X. 

In the following, we are interested in the behaviour of densities (I7j) in the ground state 
of the system. In the strong coupling limit the dominant behaviour of the ground state is 
determined by Hps- This means that the suitable fermionic coordinates are the coordinates 
which diagonalise Hps- Accordingly, we rewrite the Hamiltonian ([T]) in terms of Xj's, f//'s, 
dl's and rf-^'s defined by (jl]), expressing derivatives as [5|: 

JL = y-t,.ct.A+ y- ^Mc(rj_+ y- IHEL,^ r (9) 

fix'' Z^'^k^^ f)r. r -r ATP ^ r -r " ^ ' 

This transformation of coordinates induces a nontrivial measure in the definition of scalar 
product of the states. After performing similarity transformation and defining new Hamil- 
tonian H = SHS~^ with S being S = Yli^j ~ ^i)^^^' scalar product of new states 
|\E') = 5*1 \E') is defined with respect to desirable trivial measure. Finally, the Hamiltonian 
([1]) can be recast into the following form 



k k k m,n,m^n 



where we defined for n ^ m 



Ll = U?^-\{dld--d-dl). (11) 



The operators L"^ = U^JlJJJ" are transformed "angular momentum" generators which 
commute with Hamiltonian ([T]) and generate unitary group. Assuming that the last term in 
the Hamiltonian (fTOl) may be ignored, assumption to be justified latter on, the eigenstates 
of the Hamiltonian (fTOj) can be written as^: 

|§) = 0(a;i, a;„, yi, yN-n)d\ ■ ■ ■ 4|0)/> (12) 



^ The condition that the eigenstate of the Hamiltonian H is antisymmetric with respect to exchange of two 
fermions imphes that state transformed by aforementioned similarity transformation is symmetric with 
respect to this exchange. By this transformation the matrix degrees of freedom are effectively described 
as fermions. 
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where coordinates Xi,i = 1, . . . , n correspond to occupied, while coordinates yj = Xn+j,j = 
1, . . . , N — n correspond to unoccupied levels. In the following, we require that function 
is antisymmetric under the exchange of two indices of the occupied levels and separately 
under the exchange of two indices of the unoccupied levels. Therefore, the proper state of 
the system is obtained by antisymmetrising the state fll2p with respect to exchange of indices 
of occupied and unoccupied levels. The action of the Hamiltonian ffTOj) . without the last 
term, on the state f|T2l) reduces to the action of the Hamiltonian i^red. on the wave function 
0, with Hred. givcu as: 



k ~ k ~ k "'^f^ ~ k 



where Xi = Xi_ + g/^u"^, Vi = Vi — gj^ij^. The properly symmetrised ground state of the 
Hamiltonian -ffred. is 

~ n^^^ "^j^)^ n^^^ "^j)^"^"^^'*"^'^^'^'- (1^) 

Suppose now that we have exact ground state |\[^ exact)- Due to singularity for Xi ~ yj of the 
last term in the Hamiltonian (fTOl) . the exact ground state contains the prefactor J^- ^ [xi — yj) 
with the suitable power, as usual in the Calogero-like models j8|. In the leading order in 
g the introduction of this prefactor into the state f|T^ results in multiplication of the state 
by a constant, and formally the appearance of this prefactor is out of the scope of present 
approximation. However, the important effect of this prefactor is that the contribution to the 
expectation value of the last term in the Hamiltonian ( ITOi) in the exact ground state reduces 
to the principal value integral avoiding singularity. Using the principal value prescription we 
can expand the integrand into powers of g, showing that this contribution in the state f[T^ 
is of order g'"^ relative to the other terms in the Hamiltonian. This justifies our assumption 
that the last term of the Hamiltonian (fTO!) may be ignored and the bosonic part of the 
ground state separates into the product of the part depending solely on the occupied levels 
and the part depending on the unoccupied levels. Using the usual methods from random 
matrix theory |9|, evaluation of the densities ([7]) in approximate ground state gives 

n = e-('^^y± "'^^J^;^^'>\ (16) 

i=l 



i=l 



where Hn{x) is Hermite polynomial of order n. In the limit N — n,n oo, densities 
( IT5|) reduce to Wigner's semicircle laws, as indicated in Fig.([T]), with centres separated by 
g/uj'^ in accordance with two-cut solution previously found [3|. This result shows that in the 
strong coupling limit we still have well separated densities of occupied and unoccupied levels, 
although the probability nature of the levels dynamics allows mixing of levels. Furthermore, 
Eq. (fT5l) shows that in the case of finite number of levels and finite number of fermions there 
exist states in the previously found homolumo gap [sl]. Note that in precisely this regime, 
i.e., in the limit of strong interaction, the model defined by ([Tj) was used as a toy model for 
QCD in the analysis of the spectrum of mesons and baryons |5|. 



5 




FIG. 1: Smoothing out the densities on the upper graph is plotted p°'^'^{n = 20) and p^'^°'^'^(^N— 
n = 15), and on the lower graph is plotted /9°'^'^(n = 100) and p'^^°'^'^{N —n = S{)), with g/2uj'^ = 15. 

In the case of weak interaction, up to first order in g we can write the Hamiltonian ([1]) as 

H = e'^fHse-'^f", (16) 

where 



Lk 



The eigenstate of the Hamiltonian (1161) is 



JH 



fb 



(17) 



(18) 



where |$) is an eigenstate of Hb- The density of occupied eigenvalues in this state is 

i 

= J2{'^\'^~'^^'4d'5{x - Xi)e'"f'^\^). (19) 

i 

Expanding in g up to first order we obtain 

(20) 

where (Po'^^(x)) is the density of occupied levels for g = 0. Analogously we find 



2uj' 



(21) 



The case of weak interaction shows that starting with the system of bosons determined by 
the Hamiltonian Hb, introduction of boson-fermion interaction results in the displacement 
of the occupied levels by —g/2uj'^ and unoccupied levels by g/2u!^, whatever these densities 
are in the case g = 0. Moreover, one can show that even to the second order the boson- 
fermion interaction results only in the displacement of the of occupied levels by —g/2uj'^ and 
unoccupied levels by g/2uj'^. As this is exactly the result of the strong coupling approxima- 
tion, one might speculate that this would be also true for the intermediate g. In that case, 
in the transition region where densities, although displaced, still overlap, as in Fig([2]), the 
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FIG. 2: We sketched development of gap in the weak coupUng Hmit for N — n,n ^ oo. 

total density of levels in the large N limit (the sum of two semicircle distributions) would 
display two (ignoring the end-points) homolumo 'kinks'. 

The appearance of the gap in our picture depends on the frequency uj, the interaction 
strength g and the number of levels and number of fermions n. The control over these 
parameters effectively controls the size of the gap. A possible disappearance of the gap, 
observed in reality, can be also interpreted within our picture as a consequence of the exis- 
tence of the several subsystems whose gaps are arranged in appropriate way. Namely, one 
could introduce different frequencies Ui in bosonic part of the model, which would result in 
different gap for each subsystem, and these could be arranged to overlap. This could provide 
a simple model for exploring the insertion of a single level into the gap. 

Furthermore, the Hamiltonian ([T]) looks exactly as the Hamiltonian used to derive the 
linear Jahn- Teller effect which says that due to the degeneracy of orbital fermions levels 
the symmetry of the molecule is broken in the ground state. Although in description of 
this effect the Hamiltonian ([1]) is defined in configuration rather than in momentum/energy 
space, the homolumo gap we observe can be interpreted as manifestation of breaking of 
gauge symmetry, representing generalisation of Jahn- Teller effect. 

In our approach, the crucial point was use of suitable unitary transformation that enabled 
us to diagonalize the model in both strong and weak coupling approximation. On the 
other hand, using the conserved " angular momentum" generators and appropriate spectrum 
generating algebra, one should attempt to construct exact eigenstates of the system. We 
hope to report on the related results in future publications. 
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